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Abstract. 

We study the finite temperature free energy and fermion number for Dirac 
fields in a one-dimensional spatial segment, under local boundary conditions, corn- 
ed ■ patible with the presence of a spectral asymmetry. We discuss in detail the con- 
^ ' tribution of this part of the spectrum to the determinant. We evaluate the finite 
C ~ ) ' temperature properties of the theory for arbitrary values of the chemical poten- 
tial. 

Subject Classification 
PACS: 11.10. Wx, 02.30.Sa 

in 
o 

S3 



On 



c3 



1. Introduction 



When the Euclidean Dirac operator is considered on cvcn-dimcnsional compact 
manifolds with boundary, its domain can be determined through a family of local 
boundary conditions which define a self-adjoint boundary problem .Ij (the particular 
case of two-dimensional manifolds was first studied in [2]). The whole family 
is characterized by a real parameter 9, which can be interpreted as an analytic 
continuation of the well known parameter in gauge theories. These boundary 
conditions can be considered to be the natural counterpart in Euclidean space of 
the well known chiral bag boundary conditions. 

One salient characteristic of these local boundary conditions is the generation of 
an asymmetry in the spectrum of the Dirac operator. For the particular case of two- 
dimensional product manifolds, such asymmetry was shown, in 3_, to be determined 
by the asymmetry of the boundary spectrum. For other recent work on chiral bag 
boundary conditions, see 0H5JH>]- 

In a previous paper 0, we studied a theory of Dirac fields in one spatial 
dimension and evaluated its finite-temperature properties for two particular values of 
the parameter 9, and for restricted ranges of the chemical potential /i, since our aim 
was to (partially) answer the question posed in [H] , as to whether the fermion number 

t Member of CONICET 
§ Member of CONICET 



Dirac operator with bag boundary conditions 



2 



is modified by temperature in low dimensional bags. Here, we generalize such results 
to other values of \x. Such generalization could be of interest in the study of boundary 
effects in low-dimesional condensed matter systems and also in the treatment of open 
string theories with a nontrivial twist (analytic extension of /i) of the world sheet 

mum. 

In section |21 we determine the spectrum of the Euclidean Dirac operator at finite 
temperature for 9 = 0. With this spectrum at hand we perform, in section the 
calculation of the partition function via zeta function regularization for different ranges 
of the chemical potential. 

Section 0] is devoted to the evaluation of the free energy and the mean fermion 
number, both at finite and zero temperature. 

2. Spectrum of the Euclidean Dirac operator 

In order to study the effect of temperature, we consider a two-dimensional Euclidean 
space, with the metric (+,+). We take the Euclidean gamma matrices to be 70 = 0i, 
71 = 02 ■ Thus, the Euclidean action is 



S E = J d 2 x^{ip-A)^. (1) 

The partition function is given by 

logZ = log det(ip- A) BC . (2) 

Here, BC stands for antiperiodic boundary conditions in the "time" direction 
(0 < xo < /3, with (3 = y) and, in the "space" direction (0 < x\ < L), 



^(l + 7o)^ 



= 





^(l+7o)* 



= 0. (3) 

L 

We will follow jl lj in introducing the chemical potential as an imaginary Aq = 

—ifi. 

In order to evaluate the partition function in the zeta regularization approach, 
we first determine the eigenf unctions, and the corresponding eigenvalues, of the Dirac 
operator 

{if)+ijoH)ty = ujty . (4) 
To satisfy antiperiodic boundary conditions in the xo direction, we expand 

*(x ,x 1 )=Y / e lXX0 ^i)i (5) 

A 

with 

A ; = (2Z + 1)^, I = -00, ....oo. (6) 

After doing so, and writing ib(xi) = ( ^[ Xl \ ] we have, for each L 
(-Aj + di)x = ulp 

{-h-dx)y = u)x- (7) 
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A; = A; - i/J, = (21 + 1)~ - i^ . 



(8) 



It is easy to see that, with the boundary condition in equation @, no zero mode 
appears. For u> ^ one has, from (0), 

dfcp = —K 2 (p 

(9) 



X = --(A; + di)<p, 



where k 2 = uj 2 — Xf. 

For k 0, one has for the eigenvalues 



uj n ,i — ±y (~^J witn n=l,...,oo, / = — oo, ...,oo. (10) 

This part of the spectrum is symmetric. In the case k = one has a set of x\- 
independent eigenfunctions, corresponding to 

w,=A,. (11) 

It is to be noted that u>i — —A; are not eigenvalues. 

3. Partition function 



In order to obtain the partition function, as defined in equation @ , we must consider 
the contributions to log Z coming from both types of eigenvalues (equations (|10fl and 

CP) 



1 ds 



s=0 



where 



Ci( S ) = (i + (-ir s ) J2 

n = 1 

I = — oc 

A - — 

2 ds 



and 



\oLJ 



0»(«) 



(2i+l)^-i^ 



(12) 



(13) 



s=0 



where 



&« = E 



Z= — oo 



(14) 



(15) 



As usual, a is a parameter with dimensions of mass, introduced to render the 
zeta function dimensionless. 

The analytic extension of £2 requires a careful selection of the cut in the w-plane 
(for details, see reference [7])- The result is 



Ca(*) 



2tt 

a/3 

/3 a 



1 i/x/3 



E 



i=0 



- ( ' + 2 } -^ 



1 ifi/3 



' 2 2tt 



(16) 
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where Ch(s,x) is the Hurwitz zeta function. 

The analytic extension of £1 leads to a separate consideration of different /Lt-ranges, 
determined by the energies of the zero-temperature problem [7] . We will perform the 
extension in two of these ranges. The generalization to other ranges will become 
evident from these two cases. 



3.1. H < i 

This is the range already studied by us. For details, see reference 0. In this range, 
H13JI can be written in terms of its Mellin transform as 



(1 + (-!)-*) 

r(|) 



f v e -'[(2*) a +(c^a-*s) a ] 

Jo 



(17) 



n = 1 

/ = — oo 



This can also be written as 

(i + (-ir s ) 



&(«) = 



E 



—t-K 



e 



(v^) s r(f) — 

It ( 1 ifi 



4t 



[(^) 2 +(^-^) 2 ] 



where we have used the definition of the Jacobi theta function 

To proceed, we use the inversion formula for the Jacobi function, 
03(z, x) = ■^e^ _ ^ _ ^©3 (^r, i), and perform the integration over t, thus getting 



(18) 



Ci(«) = 



(1 + (-!)-)/? 

2ar«(0F)T(§) 



S — 1 \ 7T 2 



i=r C«(«-i)+ 



,1=1 



cosh (p0l)K s -i 



nln/3 
L 



(19) 



where £r is the Riemann zeta function. 

From (|19|l and Ijl6(l both contributions to log Z in this range of can be obtained. 
They are given by 



/?7T 

12L 



E 1 °g( 1 



e + 2cosh(fj,(3)e ^ \ 



(20) 



and 
A, = 



1 ifi/3 



2n 



1 ifi/3 



(ji 0,--f +C 5 0,- + f + insignia 0, - + 



2tt 



1 



2tt 



log 2 + log cosh 



(21) 



Putting both pieces together, we finally have 

(3tt 
12L 



„ oo 

logZ= - +^log(l + e _2z ? fi +2cosh(^/?) e - Ii £ £ ) 



n=l 

log 2 + log cosh 



|m|/3 

2 



(22) 
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3.2. f < M < X 
Again, we have 

(1 + (-!)-*) 



Ci(*) = 



r(f) 



dtt^ 



n 



I 



oc 

E *- 
1 

-co 



*[(s) 2 +(( 2 ;+i)^-^) 2 ] 



(23) 



However, in this range of fi, the contribution to the zeta function due to n = 1 
must be analytically extended in a different way. In fact, the expression cannot be 
written in terms of a unique Mellin transform, since its real part is not always positive 
(note, in connection with this that, for n = 1, eq. (|19|) diverges). Instead, it can be 
written as a product of two Mellin transforms 



dtti- 1 e~^ {2l+1) f^ +i ^ t 



dz z 2 ' 



-[(2i+l)f-V-i£]* + 



(24) 



or, after changing variables according to t' = t — z; z' = t + z, performing the integral 
over t' and the sum over I 



, n=1 , (l + (-l)- s )v^ 



2a- s r(§) 



dz ; 



fl 



(25) 



— y L / sinh(^) 

Now, the integral in this expression diverges at z — 0. In order to isolate such 
divergence, we add and subtract the first term in the series expansion of the Bessel 
function, thus getting the following two pieces 



(i + (-i)- s )V^ f*\^ 



and 



4a- s r(f + 1 
dz z~ 
fi — > —fi , 

(i + (-i)-«)V^ 



J^z) 



a lfJ,Z 



sinh(f) 



■row 2 «a-«r(f)r(2^) 7 



<iz ; 



fl 



-ft. 



(26) 



(27) 



sinh(^, 

The contribution of equation 1)26(1 to the partition function can be easily evaluated 
by noticing that the factor multiplying s is finite at s = 0. Thus, one has 

pi/12 



a n=l 
A 1.(D 



cos ( — z) 



1 



sinh(f) 



(28) 



where we have used that J_i(j^z) = y 7x7 cos (x z )- Now, in the term with fi 
one can change z — > — z to obtain 



A i,(i) - _ 



G?Z , 



COS (yZJ 



sinh(^) 



(29) 
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This last integral is easy to evaluate in the complex plane, by carefully taking into 
account the sign of fx, as well as well as the fact that j- < in closing the integration 
path, to obtain 



A"; 



n=l 



(1) 



1=1 



(30) 



or, after summing the series 

A^ 1 , = {log (l + e~ 2 ^ + 2 cosh (j-0)t 

+ | M |/3-21og^2cosh(^)H . (31) 

In order to get the contribution coming from i|27|) . the integral can be evaluated 
for SRs > 1, which gives 

(i + (-i)- s )r( s )^( a /?) s 



(27r) s 2-ir(f)r(^i) 

Ch(s, -(1 )) + Ch(s, r(l + )) 

I 7T Z 7T 



(32) 



Its contribution to the partition function can now be obtained by using that 
Ch (0, |(1 - Hf)) + Ch (0, 1(1 + = and the well known value of -£j s=0 Ch (a, x) 
|12j . to obtain 



a n=l 
A l,(2) 



21og(2cosh(^)). 



(33) 



Summing up the contributions in equations (|21|l . I|31|) and <|33[1 . as well as the 
contribution coming from n > 2, evaluated as in the previous subsection, one gets for 
the partition function 



logZ 



log! 2 cosh ( ^ 



log (l + e" 

oo 

EM 1 



2 cosh (^e" 1 ^ 



12- 1 



-2nfB 



2 cosh (/x/3)e 



-nf/3 



(34) 



n=2 



At first sight, this result looks different from the one corresponding to < j- 
(equation (|22|l ). However, it is easy to see that both expressions coincide, the only 
difference being that the zero-temperature limit is explicitly isolated from finite- 
temperature corrections. Similar calculations lead to the same conclusion for other 
ranges of fi. In those cases where /i coincides exactly with one energy level, the result 
can be shown to be the same, but series such as those in equations Ijl9|l and 1|30|) are 
only conditionally convergent. 



4. Free energy and fermion number 

From the results in the previous section, the free energy can be readily obtained. It 
is given by 
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W 



log 2 + log cosh ( — 



It is continuous, in particular, at \fi\ 
limit one has 

„,kn , (k + l)7i\ 
F(-<|,|<L_^) 

The fermion number is obtained as 
1 d \ogZ 



(35) 

= kj, k — 0, ...,00. In the low-temperature 



12L 



+ k(k + l) 



2L 




(36) 
(37) 



-H0 



sign(fx) 



(38) 



Note that it is not defined for /i = 0, where both lateral limits differ. This 
originates from the indetermination of the phase of the determinant (equation 
From a physical point of view, this reflects the fact that the sign of fi distinguishes 
particles from antiparticles. 

Particularly interesting is the discontinuous behavior of N in the zero-temperature 
limit, where y, is to be interpreted as a Fermi energy. In such limit, one has 

1W 

' -^/3^oo -ksign(fi) , (39) 



, T .fc7T (k - 



which coincides with the derivative of equation (|36|) . and is consistent with Fermi 
statistics. For /i equal to an energy level, both lateral limits differ, and -ZV(|/i| = ^) 
is undefined. 
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